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Transformation thermodynamics, as one of the important branches among the extensions of 
transformation optics, has attracted plentiful attentions and interests recently. The result of transformation 
thermodynamics, or called as "thermal cloak", can realize isothermal region and hide objects from heat. In 
this paper, we presented the concept of "reverse thermal cloak" to correspond to the thermal cloak and made 
a simple engineering definition to identify them. By full-wave simulations, we verified that the reverse 
thermal cloak can concentrate heat and realize local heating. The performance of local heating depends on 
the anisotropic dispersion of the cloaking layer's thermal conductivity. Three-dimensional finite element 
simulations demonstrated that the reverse thermal cloak can be used to heat up objects. Besides 
pre-engineered metamaterials, such reverse thermal cloak can even be realized with homogenous materials 
by alternating spoke-like structure or Hashin coated-sphere structure. 

With the born of transformation optics in 2006 1,2 , it grows prosperously with a lot of attention and has 
turned into a subject of great interest since it not only provides a general method in coordinate 
transformation but also provides an effective design tool for manipulating electromagnetic wave 3 . In 
the past few years, transformation optics has been well developed beyond the domain of electromagnetic wave to 
plentiful extensions in other fields, including acoustic wave 4 ' 5 , liquid wave 6 , matter wave 7 , and elastic wave 8 ' 9 , etc. 
The results of these transformations, or called as cloaks, can be realized by pre-engineered metamaterials or even 
homogeneous natural materials with special-designed structures. One common type of cloaks is the core- shell 
structure, which maps a concealment region into a surrounding shell region with inhomogeneous distribution of 
constitutive parameters 10 " 12 . The core-shell structure, like the typical Hashin coated sphere structure, was initia- 
tively described by Hashin 913 as a model of composite comprised of spherical grains of one phase embedded in a 
matrix of another phase. Hashin pointed out how to design the phase materials and the corresponding fractions so 
that inserting the coated sphere into the matrix would not disturb the uniform electromagnetic field outside 14 . 

Analogy to transformation optics, Guenneau et al 15 first extended the transformation methodology to the area 
of thermodynamics and proposed the concept of thermal cloak to channel the heat flow to realize a local 
isothermal region where the temperature gradient vanishes. Further studies 16 " 19 demonstrated that thermal cloak 
can be realized with homogenous and finite thermal conductivity by simulations and experiments. More spe- 
cifically, Han et al. 16 presented that thermal cloak can be realized by isotropic natural materials; Schittny et al. 17 
designed and manufactured a thermal cloak comprised of copper rings interlaid with plastic. From the view of 
guiding heat flux, Narayana and Sato 18 did excellent and systematic work about controlling the heat flow by the 
engineering materials. By shielding, concentrating, and inverting the heat current, they experimentally showed 
that artificially engineered thermal materials can yield characteristics that go far beyond the conventional media. 

The concept of thermal cloak, in fact, is derived from the optical invisible cloak, and it is aimed at hiding object 
from heat or rendering heat invisible in certain region. The thermal cloak can be achieved by anisotropic thermal 
conductivity distribution where the azimuthal component k' q is much larger than the radial component K r r 
(K' e >K' r ). Inspired by the previous work, one possible idea turns up that could we apply the thermal cloak reversely 
(K f r >K f e ) and realize a local heating region where the temperature gradient is large enough to heat up objects? 

By utilizing the transformation thermodynamics and tailoring the dispersion of thermal conductivity, here we 
demonstrated that the reverse application of thermal cloak can efficiently realize a local heating region where the 
temperature gradient is more considerable than that of the outside region. Steady full-wave simulations show that 
the temperature profile in the central region strongly depends on the anisotropy of thermal conductivity. 
Transient simulations illustrate that heat diffuses faster in the left part of the reverse thermal cloak and slower 
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Figure 1 | Temperature profile comparison among a homogenous plane, an annular thermal cloak and an annular reverse thermal cloak with constant 
boundary temperatures at 400 K and 300 K respectively: (a) a plane and the corresponding temperature curve of section line at y = 0; (b) a thermal cloak 
and the corresponding temperature curve of section line at y=0, {C=k' y =0.1, k'q = 10); (c) a reverse thermal cloak and the corresponding temperature 
curve of section line at y=0 (C=k' y =10, ic^=0.1). Isotherms are illustrated in green color in panel. Cyan dash lines are the section lines. 



in the right part, thus large temperature gradient is formed in the 
central region. By modeling a substrate consisting of reverse thermal 
cloak and assembling a cubic probe onto the central region, the 
corresponding temperature fields on the cubic probe show that local 
heating effect was enhanced with reverse thermal cloak. Two possible 
methods were also presented to realize such reverse thermal cloak. 

Results 

With transformation thermodynamics, we derived that the radial 
thermal conductivity k\ and the azimuthal thermal conductivity k'q 
are the inverse of each other, as shown in Eq. (1). 

<4=i (i) 

By defining k\ = 1 / k' q = C, where C is a constant, a cloak can 
be realized with constant thermal conductivity. For thermal 
cloaks, the constant C is usually less than l 16 , implying the radial 



conductivity k' y is much smaller than the azimuthal conductivity k'q. 
The performance characteristics of thermal cloak are totally attrib- 
uted to its anisotropic distribution of thermal conductivity. In this 
paper, we proposed the concept of "reverse thermal cloak" where k' y 
is much larger than k'q, i.e. the constant C is larger than 1. 

With full-wave simulations based on finite element method 
(FEM), the steady-state performances of the cloaks were examined. 
The large commercial software package COMSOL Multiphysics was 
used to conduct the simulations. For comparison, a homogeneous 
plane with constant thermal conductivity was also simulated. In the 
simulations, the left and right boundaries were kept at constant 
temperature, i.e. 400 K and 300 K, respectively. The top and bottom 
boundaries are thermally insulated. The interior and exterior radii of 
the annular thermal cloak (reverse thermal cloak) are 1 cm and 2 cm, 
respectively. The dimensions of the whole domain are 10X10 cm 2 . 
The thermal conductivities of the core region and the substrate are 
kept as 0.1 W/(nrK) and 1.0 W/(nrK), respectively. 
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Figure 2 | Visual illustration of the temperature profiles of (a) a thermal cloak and (b) a reverse thermal cloak. 



Fig. 1 shows the temperature profiles and corresponding temper- 
ature curves of the section lines on a homogenous plane, a thermal 
cloak and a reverse thermal cloak, respectively. From Fig. 1(a), the 
temperature profile of the homogenous plane is uniform and the 
temperature gradient is constant throughout the plane. From 
Fig. 1(b), the temperature profiles outside the thermal cloak are all 
the same with Fig. 1(a), but the heat tends to be away from the inner 
core region and travels close to the outside edge of the thermal cloak 
(C=K f r = 0.1 W/(m*K), 4 = 10 W/(m*K)). From both the temper- 
ature contour and the temperature gradient of the section line, we 
can see that thermal cloak can realize an isothermal region in the 
central domain. This implies that any object placed in the interior 
region appears to be concealed since the outside heat flow looks like 
as if there is no disturbance at all in the central region (the temper- 
ature gradient outside the thermal cloak is the same with that in 
Fig. 1(a)). This is the reason and origin of the definition of "thermal 
cloak". Fig. 1(c) shows that temperature field of the reverse thermal 
cloak in which we adjust the constant C to 10 (C= K f r = 10 W/(m* K), 
k' 6 = 0.1 W/(m*K)). We can see that the temperature profiles bend 
toward the central region and the temperature gradient in the central 
region is enhanced greatly. More vivid comparison between the ther- 
mal cloak and the reverse thermal cloak can be found in Fig. 2. The 
reverse thermal cloak tends to concentrate the heat and have the 
most heat diffuse from it, resulting in local heating and large tem- 
perature gradient. Isothermal region forms in the thermal cloak 
while a sharp temperature drop arises in the reverse thermal cloak. 
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Figure 3 | Temperature profiles of a cut plane (x= —0.5 cm) on an 
annular reverse thermal cloak with changing the constant Cfrom 5 to 20. 



To further demonstrate the local heating phenomenon, we plotted 
the temperature curves on the cut plane (x=— 0.5 cm) in Fig. 3. 
Obvious local high temperature can be observed in the central region. 
We can also see that when the constant C increases from 5 to 20, the 
local temperature in the central region is enhanced. The larger 
the constant C is, the larger the thermal conductivity anisotropy is. 
The trend in Fig. 3 implies that large anisotropy of the thermal 
conductivity distribution leads to enhanced local heating effect. 

Transient simulations of the reverse thermal cloak were carried 
out and the heat diffusion processes are shown in Fig. 4 
{C=k' t =10 W/(m*K), k' 6 = 01 W/(nrK)). In the beginning, heat 
has not conducted to the right part yet and the isotherms squeeze 
in the left part. With time elapse, heat conducts to the right part 
gradually and the isotherms stretch to the whole domain. At t=60 
second, the thermal system almost get its equilibrium state and the 
corresponding isotherms are symmetric in the left and right parts. 
No matter in the transient state (see Fig.4) or in the steady state (see 
Fig. 1 and Fig. 2), we can discover that the isotherms bend toward the 
central region in both the left and right parts. This implied that in the 
reverse thermal cloak, heat conducts faster in the left part and slower 
in the right part. As a result, compact isotherms and large temper- 
ature gradient form in the central region. 

From above analyses, we can observe the local heating effect and 
there exists a large temperature gradient in the central region of the 
reverse thermal cloak. Since the driving force of heat diffuse is the 
temperature gradient, one idea is that we could utilize the large 
temperature gradient to heat up objects. To verify our idea, we built 
a three-dimensional model in which we modeled a substrate consist- 
ing of a reverse thermal cloak (C=K f r =10 W/(m*K), k^ = 0.1 W/ 
(m*K)) and then mounted a cubic probe on the central region, as 
shown in Fig. 5(a). The dimensions of the substrate are 10X 10X0.5 
cm 3 . The dimensions of the cubic probe are 1X1X0.5 cm 3 . The 
interior and exterior radii of the annular reverse thermal cloak are 
1 cm and 2 cm, respectively. For comparison, we also built the same- 
dimensional model except that the substrate is a homogeneous plane 
without reverse thermal cloak. After steady-state simulations, we 
selected the temperature contours on the front surface of the cubic 
probe (on the plane of x=— 0.5 cm) to demonstrate our idea. The 
temperature comparisons are shown in Fig. 5(b) and Fig. 5(c). We 
can see that the maximum temperature of the cubic probe is 370 K 
when the substrate is with a reverse thermal cloak; the maximum 
temperature is 357 K when the substrate is without a reverse thermal 
cloak. Thus the reverse thermal cloak can be utilized to heat up 
objects. 

To realize such a reverse thermal cloak whose radial thermal con- 
ductivity is smaller than its azimuthal thermal conductivity, we 
designed a composite structure made by alternating spoke-like layers 
of two homogeneous materials. As shown in Fig. 6, the structure 
consists of two isotropic materials, i.e. material A and material B. 
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r =30 s f=40s t=60 s 



Figure 4 | Transient heat conduction through an annular reverse thermal cloak with the left and right boundaries at constant temperature. Snapshots 
of temperature profile with time elapse from t=5 s to t=60 s. Isotherms are illustrated with green color in panel. 
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(a) 



(b) 



Figure 6 | Simulation of a reverse thermal cloak consisting of alternating spoke-like layers of two homogeneous materials (A and B): (a) the 
temperature profile and (b) its structure. 



The thermal conductivities of these two materials are 10 W/(m*K) 
and 0.1 W/(m* K), respectively. All the boundary conditions are kept 
the same. The simulated temperature profile is similar to above 
simulations where compact isotherms and local heating effect form 
in the central region. This implies that the reverse thermal cloak not 
only can be realized by pre-engineered metamaterials, but also can be 
realized by homogeneous materials with alternating spoke-like 
structures. 

Another method to realize such kind of reverse thermal cloak with 
homogeneous materials can ascend to the Hashin coated sphere. 
According to the two-dimensional Hashin- Shtrikman formula 14 , 
the effective thermal conductivity k* can be calculated as 




Figure 7 | Simulation of another reverse thermal cloak comprised of 
homogenous materials. The thermal conductivity of the reverse thermal 
cloak is 0.918 W/(m*K), which is designed according to the Hashin- 
Shtrikman theory. The isotherms and streamlines are also illustrated with 
green and white colors in panel, respectively. 



fi = l-f2=(n/r 2 ) 2 (3) 

where K\, k 2 are the thermal conductivities of the core region and the 
shell region, respectively, fx and/ 2 are the two-dimensional volume 
fraction, Vi and r 2 are the exterior radii of the circular core and the 
cloaking shell. When k* equals to the thermal conductivity of the 
outside matrix /c 0 , reverse thermal cloak can be realized. 

According to dimensions and thermal conductivities mentioned 
above, the thermal conductivity of the cloaking shell k 2 can be cal- 
culated as 0.918 W/(m*K) in this paper. With the same method, we 
simulated the temperature field and the result was plotted in Fig. 7. It 
is seen that the green isotherms and the white streamlines are ortho- 
metric to each other in panel. The isotherms bend to the core region 
and a large temperature gradient forms there. Therefore, with the 
Hashin-Shtrikman formula, we can design the dimensions and ther- 
mal conductivities to realize reverse thermal cloaks. But when com- 
paring Fig. 7 with Fig. 4 and Fig. 6, it is also seen that the thermal 
gradient in Hashin coated sphere is not as compact as that realized by 
pre-engineered metamaterials, nor homogeneous materials with 
alternating spoke-like structures. Thus it is perceived that the local 
heating effect in Hashin coated sphere is not as strong as the other 
two, and the heat-concentration ability of Hashin coated sphere is 
also limited. 

Discussion 

In summary, we demonstrated in this study that reverse thermal 
cloak can be applied to realize local heating. Unlike the thermal cloak 
can realize isothermal region and render heat invisible from the 
outside temperature field, the reverse thermal cloak can bend the 
isotherms to the central region and realize large temperature gradient 
there. Due to the large temperature gradient, local heating phenom- 
enon was observed. Such "local heating" phenomenon can be 
enhanced by re-distributing the dispersion of thermal conductivity 
in the domain of the reverse thermal cloak. Here, a way to realize 
local heating was proposed by the reverse thermal cloak. Like the 
fabrication of thermal cloak, the reverse thermal cloak can also be 
fabricated with pre-engineered metamaterial or even with homogen- 
ous materials by spoke-like structures or Hashin coated- sphere 
structure. With the reverse thermal cloak, we can realize the local 
heating in many fields. For instance, in the electronic packaging, with 
local heating, we may avoid heating some temperature-sensitive 
components or causing printed circuit board (PCB) warping. 

It should be pointed out that the reverse thermal cloak was pro- 
posed to correspond to the thermal cloak during the research of 
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thermal cloak. In this paper, we simply defined the difference 
between these two cloaks as that: for thermal cloak, the constant C 
(= radial thermal conductivity k' y ) is less than 1; for the reverse 
thermal cloak, the constant C is larger than 1. The essence of both 
thermal cloak and reverse thermal cloak is the anisotropic dispersion 
of thermal conductivity. When compared with the thermal concen- 
trator 15,18 , the present reverse thermal cloak exhibits similar charac- 
teristics but in a much simpler engineering way. The thermal 
concentrator has specific formula to identify the dispersion of ther- 
mal conductivity 15 , while in this paper we extended the limitation of 
thermal conductivity and gave a general formula in Eq. (1). 
Moreover, we paid more attentions to the application of reverse 
thermal cloak, rather than just staying in the theoretical or laboratory 
stages. What should be emphasized is that the present reverse ther- 
mal cloak can also be applied to other shape cloaks, like square, 
rectangular, elliptic, etc. 

Methods 

When a temperature gradient exists in a body, experience has shown that there is an 
energy transfer from the high-temperature region to the low-temperature region. For 
a steady state without heat source, the control equation can be written as V(/cVT) = 0, 
where k is the thermal conductivity and T is the temperature. Similar to the invari- 
ance of Maxwell equations in the transformation optics, Guenneau etal. 15 pointed out 
that the heat conduction equations can also be written in a form-invariant manner 
under coordinate transformations. Thus, in the transformed space, the heat con- 
duction equation can be re-written as 



With 



V(kVT') = 0. 



AkA 1 



det|A| 



(4) 



(5) 



where A is the Jacobian matrix with components defined as A = d{x' ,y' ,z')/ d(x,y,z). 
The Jacobian matrix characterizes the geometrical variation in the original space 
(x,y,z) and the transformed space (x',y',z'). Eqs. (2) and (3) construct the basis of 
transformation thermodynamics and we can design and manipulate the flow of heat 
by an arbitrary coordinate transformation. 

One seminal typical design of thermal cloak was consulted to the invisible cloak 
proposed by Pendry et a/. 1,10 With the simple coordinate transformation in Eq. (4), we 
can squeeze the heat flow from a central region (r<^) in the original space (r,9,z) 
into a shell region (Ri < r<^) in the transformed space {r 1 ,9' ,z'), while the rest of 
heat flow (in the region r>^) are maintained. 



(6) 



where R 1 and R 2 are the interior and exterior radii of the annular thermal cloak. 
Substituting Eq. (6) into Eqs. (4) and (5) yields Eq. (1). 
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